Formation of collimated sound beams by three-dimensional sonic crystals 
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Focusing and propagation of sound beams is of fundamental importance in several 
branches of applied acoustics, such as tomography, acoustic microscopy and imaging or 
sonar communication. To achieve optimal focusing, and to maximize spatial quality of the 
sound beams, several mechanisms have been proposed in acoustics, like the use of acoustic 
lenses^ or the design of Gaussian beam transducers.™ Recently it has become apparent that 
the materials whose properties are modulated in space, also known as sonic crystals (SCs) in 
acoustics^ or photonic crystals in optics 4 -, can modify the spatial dispersion of propagating 
waves. This feature opens new possibilities to control the diffractive broadening of sound 
beams. In particular the beams can propagate in modulated material without diffraction 
(the effect also referred to as self-collimation), as predicted and demonstrated in optics^ and 
in acoustics.- Self-collimation is based on the existence of flat segments of spatial dispersion 
curves (the curves of constant frequency in fc-space). More recently the three-dimensional 
(3D) self-collimation by SCs was experimentally demonstrated,- which is based on the 
formation of flat of the isofrequency surfaces. 

In addition to non- diffractive propagation inside the SCs, the modification of the spa- 
tial dispersion can also produce phenomena outside the crystal such as lensing 8 ' 9 and 
superlensing. 10 These beam propagation effects behind the SCs are related with the nega- 
tive diffraction inside the periodic structure. The character of the beam propagation behind 
the SC depends on the wave front of the beam acquired in the system. In particular, if 
the wave front of the beam acquires positive curvature (due to propagation in a material 
with negative, or anomalous diffraction), the beam can be focused behind the modulated 
medium, which enables above discussed lensing and superlensing effects. 

Although the focusing of sound beams behind a 2D SCs is being intensively investigated 11 
the overall picture of the beam formation and propagation is still unclear. Apart from the 
above mentioned phase transformation effects due to the negative diffraction of waves prop- 
agating inside of the SC, spatial (or angular) filtering effects also come into play. The nega- 
tively curved segments of dispersion lines are generally surrounded by the angular bandgaps, 
which are angular areas where sound cannot propagate. The latter results in a modifica- 
tion of the angular spectrum of the beams,— recently demonstrated in both optics^ and 
acoustics. 14 These two beam formation mechanisms combine, and give rich possibilities of 
formation of the beams with desired spatial characteristics (angular distributions) and with 
desired character of focusing. 
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FIG. 1. (Color online) (a) Experimental configuration, (b) unit cell of the 3D SC, (c) isofrequency 
surfaces and (d) cross-section of dispersion surfaces. 

In the present work we study, experimentally and theoretically, the sound beam formation 
behind a 3D SC with a woodpile-like structure. We experimentally demonstrate the forma- 
tion of high spatial quality and well-collimated beams behind the SC, based on the above 
described spatial filtering and negative diffraction effects. In this letter, we present first the 
isofrequency contours of SC, and identify the ranges of the frequencies where the curvature 
of the isofrequency surfaces is positive. Based on these results we design the samples, and 
perform the sound beam propagation experiments. The most important result is the exper- 
imental demonstration of the formation of a well focused beam. Finally, we investigate the 
beam propagation in a simplified paraxial approximation, and obtained a good quantitative 
interpretation of the experimental measurements. 

Figure [T^, shows the experimental setup. The 3D SC is formed by two 2D structures of 
square symmetry, embedded one into another after a relative rotation by 90 degrees, which 
results in a 3D woodpile-like structure.- Each of 2D structures consist of 20x20 matrix of 
steel cylinders of a radius r =0.8 mm, and the lattice constant a = 5.25 mm (see unit cell in 
Fig. Hb). The beam, emitted by an ultrasonic source, propagates through the SC along the 
z direction. The acoustic field is measured by a needle hydrophone positioned by a three 
motorized axes governed by acquisition system. As shown in Fig. [1^, the experimental set 
up is immersed in a plexiglass tank filled with distilled water. 

The eigenfrequency analysis of the sound wave propagation was performed numerically 
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FIG. 2. (Color online) Experimental measurements and simulations of the acoustic field pattern 
behind the SC at (a) 235 kHz, (b) 250 kHz and (c) 260 kHz. (I) XZ-cross-section of the beam 
behind the SC; experimental (2) and numerical (3) ultrasound field distribution in the XY plane 
at point zi=195 mm behind the SC. 

using Finite Element Method.— The periodicity of the system is considered by imposing 
Bloch-Floquet boundary conditions of the unit cell (Fig. db). The path around the first 
irreducible Brillouin zone represents the main directions of symmetry in 3D. We analyze the 
propagation along TX direction in the present work. 

Figure []J shows the isofrequency surfaces for three different frequencies (230, 240 and 
250 kHz) in the second band as well as the cross sections of the isofrequency surfaces by 
k z = 0, k x = and k y = planes respectively. Fig. [Tfc shows the quarter of the isofrequency 
"bubble" for these three frequencies. The isofrequency lines in k z =0, plane are shown in 
detail in Fig. [T]i. The lowest of highlighted frequencies (230 kHz) corresponds to non- 
diffractive propagation inside the SC (flat isofrequency line). The isofrequency surfaces 
(and the lines in the cross plane) at slightly higher frequencies have areas with a positive 
curvature, which cause the desired focusing behaviour. 

The experimental measurements of the beams propagating behind the SC are summarized 
in Fig. [2j Three different frequencies are represented in (a) upper, (b) middle, and (c) 
bottom panels. The upper panel (235 KHz) shows the beam propagation for frequency 
corresponding to self-collimation inside the crystal.- The bottom panels in Fig. [2] (260 
kHz) show the beam propagation for the case when a strongly curved and relatively small 
"bubble" of isofrequency surface occurs (Fig. [Tfc). As the area of the isofrequency surface 
responsible for the negative diffraction and eventually for focusing is very small, just the 
central (paraxial) part of the angular spectrum is focalized. One part of the remaining 
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FIG. 3. (Color online) (a) Density plot representing the measurements of the on-axis intensity 
distribution on the frequency-z plane (distances measured from the transducer). Continuous green 
line corresponds to analytical fit (see text). Black area represents the space occupied by the 
SC. The white dashed line represents the point Z\ (b) Continuous line represents the measured 
beam amplitude in normalized units at point Z\ and red dotted line represents the measured beam 
width, (c) Experimental profiles in the x-axis of the beam at Z\ point and (d) measured amplitude 
in normalized units on the z-axis at 250 kHz. 



angular components is reflected, as it corresponds to the angular bandgaps. The other part 
of angular components propagates along different directions, giving rise to side-lobes as seen 
in Fig. [2fc. The intermediate situation, corresponding to the frequency 250 kHz, is shown 
in the middle panels of Fig. The diffraction is negative for the propagation inside the SC 
due to positive curvature of the dispersion curves (see Figs. [Tfc and [Hi). The isofrequency 
"bubble" is large enough to transmit a sensible portion of the angular spectrum. This case 
is most relevant for the goals of this work. 

We interpret the focusing of the beam in the terms of Ref. 
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the sound beam propagat- 



ing in bulk of SC with negative diffraction accumulates the increasing positive (anomalous) 
curvature of the beam wave front. Behind the SC, the propagation in the (normally diffract- 
ing) homogeneous medium compensates the accumulated negative diffraction acquired inside 
the SC. The beam is focalized at some distance Zf where the negative diffraction inside the 
SC and positive diffraction behind the SC compensate one another. 

The analytical estimation of the focal distance is possible considering the approximation 
of small filling fraction of the SC, / = V s /V uc (V s and V uc are the volume occupied by the 
scatterer and the unit cell respectively). In this approximation the diffraction coefficient (i.e. 
the curvature of the spatial dispersion curve, and/or surface) can be analytically calculated.— 



Following the above interpretation, the negative diffraction of the SC is compensated at a 
distance Zf behind the SC (measured from the input plane of the SC): 

where L is the length of the SC, and AQ = (Q g — Q)/Q g with Q g = oj g a/2-KCh being the 
normalized Bragg frequency and Ch the speed of sound of the host medium, i.e. in water. 

The experimental results were compared with the analytical study of beam focusing, as 
shown in Fig. [3J In Fig. [3k the absolute value of the intensity behind the crystal on the z- 
axis is mapped depending on the frequency. Green continuous line represents the analytical 
fit of the focal distance calculated from Eq. ([1]) considering L = 20a, Q g = 1, and / = 0.05. 
The parameter / is a fit parameter. We notice that due to the fact that Eq. [I] has been 
obtained for 2D structures with low filling fraction, we use a fit parameter to take this into 
account. For this case, the frequency of zero diffraction point or self collimation corresponds 
to Qzdp = (1 - / 2/3 )0,, (vzdp = 238 kHz). 

The focusing for a frequency range around the optimal one is evidenced in FigsJ3b-d. 
FigfSb shows both the amplitude at point z% (black continuous line) and the beam width 
(red dots) depending on the frequency. Figs. [3fc and(3U show the profile at Z\ along the x-axis 
and the transversal cross-section along the z direction in x = for the focusing frequency 
250 kHz respectively. 

In summary, we have experimentally demonstrated the collimation of the beams behind 
a 3D SC. The obtained results are interpreted and analyzed in terms of curvatures of spatial 
dispersion curves and surfaces of the SC, and rely on the negative diffraction close to the edge 
of the propagating band. The experimental results fit well with the numerical simulations as 



well as with analytical predictions in Ref. 



The tunability of the focal distance has been 



also demonstrated, showing that the beam intensity in the focus as well as the broadening 
of the beam along the propagation depends on the frequency, which give additional options 
for applications. 

The overall focusing process is interpreted in terms of the interplay between two related 
but different effects: the focusing of the beam, due to curvature of spatial dispersion curves; 
and the spatial filtering effect, due to the size of the isofrequency "bubble". The optimum 
result comes from a compromise between these two ingredients. 
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